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Abstract 

This paper studies the performance of Mobile Ad hoc Networks (MANETs) when the nodes, that form a 
Poisson point process, selfishly choose their Medium Access Probability (MAP). We consider goodput and delay 
as the performance metric that each node is interested in optimizing taking into account the transmission energy 
costs. We introduce a pricing scheme based on the transmission energy requirements and compute the symmetric 
Nash equilibria of the game in closed form. It is shown that by appropriately pricing the nodes, the selfish behavior 
of the nodes can be used to achieve the social optimum at equilibrium. The Price of Anarchy is then analyzed for 
these games. For the game with delay based utility, we bound the price of anarchy and study the effect of the price 
factor. For the game with goodput based utility, it is shown that price of anarchy is infinite at the price factor that 
achieves the global optima. 

Index Terms 

Game Theory; Mobile Ad hoc Networks (MANETs); Pricing; Medium Access Control; Stochastic Geometry; 
Replicator Dynamics 

I. INTRODUCTION 

In this paper we study competition for network resources at the medium access control (MAC) layer. 
It is well known that computing the Nash equilibria in games is in general a hard problem. Indeed, 
this problem falls into a class of problems introduced by Christos Papadimitriou in 1994, called PPAD 
(Polynomial Parity Arguments on Directed graphs). In view of this complexity, it becomes attractive to 
identify classes of games for which one can compute the equilibria at a low complexity. We thus study 
a MAC game under some statistical assumptions on the location of the nodes, which are on one hand 
reasonable in many real scenarios, and on the other hand, allow for tractable and in several cases, even 
explicit expressions for the Nash equilibria. 
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We consider slotted time, and assume that the transmitters are synchronized. The basic assumptions on 
our model are 

• The location of the transmitters at each time slot forms a homogeneous Poisson point process (P.p.p); 

• Medium access is controlled using Aloha; 

• Transmission success is based on signal to interference and noise ratio (SINR) being larger than some 
threshold; 

• We assume saturated sources, i.e., every mobile has always a packet to send. 
The Poisson assumption means that 

• The number of mobiles in disjoint sets are independent; 

• The number of mobiles in any given set follows a Poisson distribution. This class of point processes 
maximizes entropy. It is often used for modeling the location of users in e.g. mobile ad hoc networks. 

We first consider the problem in which each node is a player: it chooses a MAP so as to maximize 
its own performance metric. We consider utility functions that model the tradeoff between quality of 
service indicators (such as goodput expected delay) and power consumption related to the transmission. 
Our Aloha assumption on the MAC protocol leads to explicit expressions for the performance metrics 
of interest, which allow us to derive several interesting properties of the network. We leave the study of 
other, less tractable, CSMA type of MAC protocols for future study. 

Our goals are to obtain symmetric Nash equilibria (SNE), possibly mixed, and to study their properties 
Q. By restricting to symmetric games, and by restricting to choices of MAP that are not functions of 
the locations of players or interferers, we have to consider the utility evaluated only at points where the 
strategies of all other players are the same, say p. We search for a p such that if all use it then a player 
who deviates unilaterally will not improve its utility. Thus utility of a player can be viewed as the utility 
of a player in a much simpler two person games, and any equilibria in the two player game is a symmetric 
equilibria in the original game. 

We are interested in particular in the case where the power consumption disutility represents a pricing 
decision of the network that wishes to determine a pricing that will induce an efficient equilibrium (in 
terms of the achieved goodput). Alternatively, the pricing may be taken such as to maximize the network's 
revenues. 

Our main findings are: 

1) Considering the goodput as the quality of service to be maximized, we observe the tragedy of the 

'We shall compute the SNE but do not address the question of how equilibrium is reached 
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commons HI: the utility at equilibrium is zero for large values of the pricing parameter. Thus the 
price of anarchy is infinite. We show however that there exists a pricing parameter for which the 
goodput at equilibrium equals the one obtained under global cooperative throughput maximization. 
We show that when each node uses the replicator dynamics to update its MAP, the system converges 
to the SNE. 

2) Considering the delay as the quality of service to be minimized, we observe that the price of anarchy 
is bounded for any price value. Here too, there exists a pricing that results in global optimum at 
equilibrium. We show that the SNE is not unique. The range of price parameters for which two 
SNE exist is characterized. 

As the price of anarchy is infinite for certain price factors, it may seem therefore that high prices have 
a negative effect on the network performance. However if we consider instead the network performance 
or the monetary profit of the operator, we discover that pricing can induce an equilibrium for which these 
measures coincide with the global optimal values. 

The Poisson assumption on the location of nodes allows us to obtain utilities in a surprisingly simple 
explicit form, which in turn allows us to obtain much insight on the property of the equilibria and on the 
role of the pricing. Such Poisson assumptions are often justified. For instance, in 10, in the context of 
cellular networks, the authors obtain explicit expressions for the coverage and throughput with the Poisson 
assumption on base station locations. Their numerical investigation show that these expressions closely 
capture the real behavior of the network. 

Related Work: There are several papers that model the nodes in the Aloha system as selfish users. In 
(31 a game theoretic model is proposed to analyze the performance of the Aloha network. It is assumed 
that when a collision occurs, some of the transmissions can be successful. The authors in [4] study the 
performance of slotted Aloha with n nodes that are selfish. Any collision results in the loss of the packet 
and there is a cost associated with unsuccessful transmissions. They characterize the equilibrium point as 
a function of n and study its limiting behavior. In [5] slotted Aloha is considered with n selfish users 
and the distributed choice of the retransmission probabilities is analyzed. It is shown numerically, that by 
adding a retransmission cost, the throughput at equilibrium equals the optimal team throughput. Pricing 
is used in the context of the power control games in [|6l to improve the performance at equilibrium. 

All these paper do not take the geometric aspects of the node location into account , which is important 
in the wireless context. We consider in this paper the geometric aspects of the MANET and analyze the 
network performance at equilibrium. Also, unlike in all other papers, where simultaneous transmissions 
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lead to transmission failure, we assume that the probability of successful transmission depends on SINR 
at the receiver. For games that take geometric aspects into consideration see 0. 

The paper is organized as follow. In Section [TT] the Poisson bipolar MANET model is set up and 
performance metrics of interest are discussed. Section [III] considers the team case in which all nodes use 
the same MAP. Section ITVl considers a scenario in which all nodes are selfish. A medium access game is 
defined among the nodes. Section |V] studies this game with a utility function that is based on the goodput 
and a pricing based on the transmission energy costs. Section [VI] studies the medium access game with 
another utility that is based on potential delay as performance metric and transmission costs. Section [VTll 
analyzes the price of anarchy for these games. We end with some remarks and future work in Section 
IVIII1 All the proofs are given in Appendices. 

II. MODEL AND SETUP 
Consider the simplified mobile ad hoc network (MANET) model called the Poisson bipolar model 
proposed in (S). Assume that each node follows the slotted version of the Aloha medium access control 
(MAC) protocol. Each dipole of the MANET consists of a transmitter and an associated receiver. Let 
the sequence {Xj,t/j}j>i denote the location of the transmitters and receivers, where y i: is the location of 
receiver associated with the transmitter at X;. We assume that the transmitters $ = {Xj}j>i are scattered 
in the Euclidian plane according to an homogeneous Pp.p of intensity A. In this paper we consider a 
scenario in which each receiver is at a fixed distance r > from its transmitter, i.e., \Xi — yi\ = r for all 
i |^|. Consider a snapshot of the location of the transmitters. Let the sequence n = 0, 1, 2, • • • denote the 
common sequence of time slots with respect to which all nodes are synchronized. We associate with each 
node a multi dimensional mark that carries information about the MAC status and the fading condition at 
each time slot. We follow the notation of [9][Chap. 17]. Let the sequence Mi = {e«(n), Fi(n)} n > denote 
the marks associated with node i, where 

• 6i = {ej(n)} n > denotes the sequence of MAC status of node i. e^n) is an indicator function that 
takes value 1 if node i decides to transmit in time slot n; otherwise it takes value zero. The random 
variable e»(n) are assumed to be i.i.d in % and n, and independent of everything else. 

• Fi — {F- (n) : j > l} n >o denotes the sequence of channel conditions between the transmitter of 
node i and all receivers (including its own receiver). It is assumed that channel conditions are i.i.d 
across the nodes and time slots, with a generic distribution denoted by F with mean 1/fi. 

• The marks are assumed to be independent in space and time. 

2 Our analysis continues to hold when the distance between transmitters and their receivers are i.i.d. 



The probability that the ith node transmits in time slot n is p : = Prje^ra) = 1} = E[ej(n)] (Medium 
Access Probability (MAP)). This defines a pair of independent Poisson processes at each time slot n , one 
representing transmitters $ 1 (n) = {Xi, ej(n) = 1} and the other non-transmitters = {JQ, ej(n) =0} 

with intensities pA and (1 — p)A respectively. All the nodes transmit at a fixed power denoted by P. 

Let l(x,y) denote the attenuation function between any two given points x,y £ M 2 . We assume that 
this function just depends on the distance between the points, i.e., \x — y\. With a slight abuse of notation 
we denote this function as l(x, y) = l(\x — y\). We assume the following form for this attenuation function 

l{x, y) = (A\x - y\)-P for A > and p > 2. (1) 

A signal transmitted by a transmitter located at Xi is successfully received in time slot n if the SINR, 
at the receiver, is larger than some threshold T, i.e., 

PFH(r) 

SINR^n) := ^ > T, (2) 

where W(n) denotes the thermal noise power at the receiver and I$ir n \ denotes the shot noise of the Pp.p 
$ 1 (n) in time slot n, namely, If,i( n )(jt/i) = Ylx-e&fa) PFl(\Xj — yi\). We assume that the noise is an i.i.d 
process. 

Consider a typical node at the origin, X = with mark M (0) = (eo(0), -Fo(O)) at n = 0. The typical 
node is said to be covered in slot n = if © holds given that it is a transmitter. Then the coverage 
probability of the typical node is 

PFl(r) 



> T 



e (0) = U , (3) 



h Ho) (y ) + W(0) 

where P° denotes the Palm distribution ifTOl rChap. I] of the stationary marked Pp.p $. Note that due 
to time-homogeneity this conditional probability does not depend on n. By using Slivnyak's theorem 
ifTOl . the coverage probability of a typical nodes when all other nodes use the same MAP is evaluated in 
[fTTTl . Continuing the notation used in [fT2l we denote this coverage probability (non-outage probability) 
as p c := p c (r,pX,T). Consider a tagged node that uses MAP p'. Then the tagged node is a transmitter 
with probability p' and a non transmitter with probability (1 — p'). We refer to the product of the MAP 
and the coverage probability of the tagged node as goodput and denote it as g(p',p) := p'p c (r,p\, T). We 
shall be interested in performance metrics of the form 

g(p\p) 1 - a 

1 — a ' 

where a > denotes the fairness parameter. In this paper we consider the performance metric 
corresponding to a = and a = 2. For a = 0, the performance metric corresponds to the goodput, 
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and for a = 2, it corresponds to the negative of potential delay introduced in IfTBI , which is defined as the 
reciprocal of the rate. Potential delay may be interpreted as the delay incurred by a node in successfully 
delivering its packets at the receiver. Indeed, it has been shown in [[9]| [Sec. 17] that the local delay is given 
as reciprocal of goodput when the node locations form an i.i.d process across time slots. In all other cases 
it is a lower bound. Local delay is defined as the number of time slots needed by a node to successfully 
transmit a packet (with retransmissions). We denote the potential delay as t(p',p) = l/g(j>',p). 

In the following sections we consider two scenarios. First we assume that the nodes of the MANET 
cooperate, i.e., use the same MAP that is assigned to them in each time slot. The value of a MAP that 
optimizes the spatial network performance is evaluated. We then consider a game scenario in which each 
node is selfish and chooses a MAP that optimizes its own performance taking into account the transmission 
costs. We study the effect of the transmission costs on the network performance at equilibrium and look 
for a cost factor that results in improved spatial network performance at equilibrium. 

III. RATE CONTROL: THE TEAM CASE 
In this section we assume that all the nodes belong to a single operator, and transmit at the MAP set 

by the operator. The following proposition immediately follows from [8J[Lemma 3.2]: 

Proposition 1: Let each node in the the Poisson bipolar transmit with MAP p and F be Rayleigh 

distributed with mean \ then the goodput is 

9(P) = pexp{- 27 rA^ _JL_ d «} x ^T/Pl(r)), 

where ip w (-) denotes the Laplace transform of the noise power W. 

Corollary 1: For F as above, zero noise power W = 0, and the path loss model in (Q~|), the goodput of 
a typical node is 

g(p) =pexp{-2n\pr 2 T 2/l3 K((3)}, (4) 
where K(J3) = mm^zlM and r(x ) = J- ^ e -». 

Hence forth we adopt the assumptions of Corollary Q] in all the subsequent calculations. However, our 
results hold when W has any distribution with finite mean, as it appears as a constant multiplicative factor 
in ©. For notational convenience we write C := C(j3,T,r) = 2i\r 2 T 2 ^ K(j3). 

The operator is interested in optimizing the social performance of the network. In particular, we assume 
that the operator aims at maximizing the mean goodput per unit area or minimizing the mean delay per unit 
area. The performance seen by a typical node can be used to derive the spatial performance of the Poisson 
bipolar MANET. Campbell's formula ifTOli rSec. 2.1.2] for stationary Poisson point processes ensures that 
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the performance experienced by a typical node is also that of the spatial average performance of the 
Poisson MANET. The mean goodput per unit area is then the product of the intensity of the P.p.p and the 
goodput, i.e., Xg(p). This quantity is referred to as the density of success and denoted by d suc (r,pX,T) 
in [9] [Chap. 16]. We denote this term simply as d s (p). Similarly, the mean delay per unit surface area is 
given by Xt(p). We denote this spatial performance metric as d t (p) and refer to it as spatial delay density. 
Note that the spatial delay density is the reciprocal of the density of success multiplied by a factor A 2 . 
Hence the MAP that maximizes the d s (p) also minimizes d t (p). The MAP that optimizes the density of 
success is given in |Hl[Prop 16.8] and [|9]|[Corol. 16.9] : 

Proposition 2: Under the assumption of Corollary [T] the MAP that maximizes the density of success 
and minimizes the density of delay is given by 

p m = min{l,l/AC}, (5) 

and the corresponding optimal density of success is given by 

f l/{eXC), if AC > 1 
d s {p m ) = < (6) 
[ Aexp{-AC}, if AC < 1, 

and the corresponding optimal delay density is given by 

A 2 eC, if AC > 1 

dt(Pm) = { (7) 
Aexp{AC}, if AC<1. 



IV. RATE CONTROL: THE NON-COOPERATIVE CASE 

In this section we assume that each node of the Poisson bipolar MANET is a selfish player. We consider 
a non-cooperative case, and model it as a game with an infinite number of players as follows. We use 
node and player interchangeably. 

Each node can take two actions: transmit (T) or no-transmit (NT). We assume that when taking the 
decision whether to transmit or not, a player does not know the positions of other mobiles nor the level 
of the SINR. A mixed strategy chosen by a node corresponds to its MAP. A player chooses its MAP once 
for all and uses always the same MAP. The choice is done without knowledge of the realization of the 
interference. 

We shall restrict to symmetric Nash equilibria (SNE) , in which all nodes use the same MAP p at 
equilibrium. The utility of a tagged player can then be written as a function of his strategy p' and of the 
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strategy p of all the others. We denote the utility of the tagged player as [/(;/, pj^. The objective of each 
player is to maximize its utility. Let U(p*,p*) denote the value of the utility function at SNE p*. 

Definition 1: p* E [0 1] is said to be a symmetric Nash equilibrium if for any node the following holds 

U(p*,p*) = max U(p,p*). 

P6[0 1] 

In the rest of this section we look for the appropriate utility functions that characterize the performance 
of an individual player. We begin with computing the goodput of a tagged node. The next lemma follows 
from Lemma 3.1 and Corollary 3.2 in (8). 

Lemma 1: Consider a Poisson MANET with the assumptions in Corollary [TJ Let a tagged node transmit 
with MAP p', while all other players transmit with MAP p. Then the goodput of the tagged node is 

g(p',p) = p' exp{-p\C}. (8) 

where C = 27rr 2 T 2 ^ K(/3) is defined earlier. 
We note the following: 

1) The goodput is monotone increasing in p' . Hence, if the objective of each mobile is to maximize 
its goodput, then the only equilibria is p — 1 for all nodes. 

2) Under the conditions of the above lemma, the expected transmission energy consumption of the 
tagged node per time slot is proportional to p' and does not depend on the MAP p of the other 
nodes. Thus the ratio between the goodput and expected energy consumption does not depend on 
p' . We conclude that any p is an equilibria when the criterion of each node is that of minimizing 
the ratio between goodput and energy. 

Keeping the above remarks in view and noting that the utility related to goodput is not a linear function 
of the MAP, see for example [fT4ll . we shall be interested in utilities] of the form 

U(p',p) = f(p',p) - w(p) (9) 

for each player, where / gives the performance measure of interest and w is related to dis-utility. We 
consider a dis-utility based on the expected transmission energy costs incurred by each node. Let p denote 
the price per unit transmission energy for each node. Then the expected transmission cost for the tagged 
node that uses MAP p' is pp'P. We define w(p') = pp\^ for a node that uses MAP p'. 

3 The game has infinitely many players, hence the utility should be defined on a infinite product strategy space. However, since we are 

restricting to the case where all players other than tagged node use the same strategy p, we write utility as a function of two arguments 
4 With some abuse of notation we use the notation U (p) to denote the expected system utility when the same p is used by all players. We 

consider utilities which attain optimal values for some p and denote the latter, again with some abuse of notation, as p m . It will be clear 

from the context if p m is an optimizer of the spatial performance metric or the team utility. 
5 without loss of generality we assume P = 1 
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Assume that the function / is a concave function in p' and continuous in p; then the arguments in 
031 [Thm. 1] can be used to show the existence of SNE. We state this result in the following lemma. The 
proof is given in Appendix IVIIIl 

Lemma 2: Assume that the utility function U(p',p) is concave in p' and continuous in p. Then SNE 
exist. 

In this paper we assume that p is set by a central agent who can regulate its value. The aim of the 
central agent is to optimize the performance of the network at equilibrium. We can also assume that the 
aim of the central agent is to maximize the network revenues at equilibrium. We refer to p as the price 
factor. 

In the next two sections we consider two utility functions defined in term of goodput and potential 
delay as the performance measure. We evaluate the MAP at equilibrium and the corresponding system 
performance. This system performance is then optimized by searching for the best price factor. The best 
achievable system performance, at equilibrium, is then compared with that evaluated when nodes act as 
a team. 

V. GOODPUT BASED UTILITY 
Assume that each node of the MANET is interested in maximizing its goodput taking into account the 

energy cost it incurs. Then by taking f(p',p) = g(p',p) in ©, we define the utility as 

U(p',p) = g(p',p) - pp' = p' {exp{-p\C} - p} . (10) 
The objective of a each node is to choose a MAP that maximizes its utility, i.e., 

p' e argmax^,^ U(p',p). 

This utility function is a linear function in p' and convex in p. Then from Lemma [2l a SNE exists. We 
proceed to calculate the SNE of this game. 

When p > 1 the slope of the utility of the tagged node is non positive irrespective of the MAP of the 
other nodes. Then the optimal strategy for each node is to choose p = 0, which is a dominant strategy 
and hence the SNE. When p < 1 consider the following two cases. 

Assume p > exp{— AC}: In this case the slope of the tagged node is always positive. Then the optimal 
strategy for the tagged node is to choose p = 1 irrespective of the MAP chosen by the other nodes. Thus 
p = 1 is a dominant strategy and hence also is the SNE. 

Assume p < exp{— AC}: If each node other than the tagged node chooses a MAP such that 

exp{- P AC}=p, (11) 



then the utility of the tagged node in (fTOb is set to zero and is not affected by its strategy, i.e., the taged 
node becomes indifferent to its own strategy. Further, if any of the nodes deviates from the MAP that 
satisfies CCD), it will not gain anything given that all other nodes set their MAP value as in (ITT]) . Hence 
the MAP satisfying (fTTT) constitutes the SNE. We summarize the above observations in the following 
proposition. 

Proposition 3: For any given C, A, and p > 

. if p > 1 then p* = is the SNE; 

. if exp{-AC} > p then p* = 1 is the SNE; 

. if exp{-AC} < p < 1 then p* = =gp is the SNE. 

The goodput of each node at equilibrium is given by 

0, ifp>l 

9(P*,P*)= I exp{-\C}, ifp<exp{-AC} (12) 

^ if exp{-AC} <p< 1. 

With the expression for goodput at equilibrium, we can look for the value of p that maximizes it. We 

assume that the objective of the central agent is to maximize the density of success at equilibrium. Then 

the optimization problem of the central agent is given by 

. . — Aplogp ,. , . , 

maximize — — subject to exp{— AC} < p < 1. (13) 

p AC 

The objective function in (fT3l) is a concave function of p > attaining its maximum at p = 1/e. If 
AC > 1 the optimal price factor lies within the constraint set and the operator can set p* = 1/e. Suppose 
AC < 1; then the objective function is decreasing in the interval exp{— AC} < p < 1 and the maximum 
is attained at p* = exp{— AC}. The maximum density of success at equilibrium with the optimal price 
factor is 

f l/(eAC), if AC > 1 
d s (p*,P*) = < (14) 
[ Aexp{-AC}, if AC < 1. 

Comparing (fl4)) with the global optimal density of successful transmissions given in ©, we have the 
following result. 

Theorem 1: The global optimal density of successful transmissions can be attained at equilibrium by 
setting the pricing factor p as follows: 

' lie, if AC > 1 

P* = { (15) 
exp{-AC}, if AC < 1. 
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Also, note that with the pricing factor p* in (fT5T ), the MAP of each node at equilibrium is the same as 
that achieving the global optimum given in ©. Thus by pricing appropriately, one can use the selfish 
behavior of the players to reach an equilibrium at which the global optimal performance is attained. 
Replicator Dynamics: In this subsection we briefly discuss how we can use tools from the population 
games and evolutionary dynamics to study aggregate behavior in the network. We can think the set of 
nodes in the Poisson bipolar MANET as a single population, where each node can take two actions. 
Let each node decide to transmit with probability p, then each node will be in state of transmission 
with probability p and in the state of non-transmission with probability probability 1-p. Then, by our 
assumption that nodes decisions are i.i.d in time and space, a fraction p of the population will be in the 
state of transmission and the other (1-p) fraction in the non-transmission state. In short we denote the 
state of the network as p when every node transmits with probability p. 

One of the most frequently used dynamics to describe the evolution of behavior in population games 
is the replicator dynamics ifToll . It describes the evolution of the fraction of the population in each state. 
In terms of population games, we can interpret the probability of successful transmission with the mean 
transmission energy costs as the fitness function for each node. Note that in our medium access game the 
fitness of a node that is not transmitting is zero. Let us denote the utility /fitness of a node that chooses 
action T while the rest of the population is in state p' as U(T,p'). The fraction of nodes that uses action 
T evolve according to the replicator dynamics as follows IfTBTi rCh. 5,6]: 

dp(t)/dt=p(t)[U(T,p(t))-U(t)} 

where U(t) denotes the mean fitness of the population. Taking (flOl) as the fitness function, the replicator 
dynamics is given by 

^ = p(t) {e-^M _ p _ p ( t)e -ACW*) _ p{t)p} = p{t){1 _ p{t)) {e -AC P ( t) _ p} 

From the above equation the stationary point is clearly or 1 or depending on the value of p, which 
is in agreement with Proposition [31 Figure |4] shows the convergence of replicator dynamics to the SNE 
starting from any interior point. 

VI. DELAY BASED UTILITY 

In this section we define the utility of each player in terms of the delay associated in delivering the 
packet successfully at its receiver, and the transmission costs. As earlier, let a tagged node incur a cost 
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of p units per unit energy dissipated. By taking f(p',p) = —t(p',p) in ®, we get the following utility 
for the tagged node: 

^''' VmUi '^ <16) 

where p' is the MAP of the tagged node, and p is the MAP used by all other nodes. 

The objective of each node is to choose a MAP that maximizes its utility function (fT6l) . Or equivalently 

it can be given by the following minimization problem: 

exp-fpAC} . , . . 

minimize h p p subject to < p < 1. 

p' p 

For a given value of p, the utility function in (fT6l) is a concave function in p' and continuous in p. Then 
by Lemma |2] SNE exist. We proceed to calculate the SNE by a direct computation. 
Differentiating the utility function with respect to p', equating to zero, and simplifying, we get p' = 
exp{p\C} / y/p, where C = C/2. This equation gives the best response of the tagged node when all other 
nodes use MAP p. If there exists a p £ [0 1] such that p' = p, then p is the SNE of the game. Hence we 
look for the conditions when the following fixed point relation holds 

expfrACfr 

VP 

Lemma 3: If XCe > ^/p, p* = 1 is the unique SNE. If XCe < VP> the SNE is 

p* = minjiw(-— j (18) 



xc V VP J 

where e denotes the base of the natural logarithm and W(-) denotes the Lambert function IfTTTl . 

Proof: By using the relation e x > xe for all x > 0, it is easy to see that p* = 1 is the SNE when 
XCe > VP- Assume XCe > yfp; rearranging the terms in (flTT) . we can write —pXC exp{— pXC} = = ^-- 
Using the definition of the Lambert function [fTTl . we obtain 



p = ^W(~J. (19) 



1 



Let us briefly recall the properties of the Lambert function in the footnote below With this explicit 
expression for the SNE we can characterize the effect of the price factor on the SNE. For some values of 

6 

• The Lambert function is two- valued in the interval [— 1/e, 0]. The two branches of the Lambert function denoted as Wo(-) and 
W-i(-) meet at —1/e and the take value —1 at this point. 

• In the interval [—1/e, 0] Wo(-) is a continuous and increasing function taking value in [—1, -co]. 

• In the interval [—1/e, 0] W_i(-) is a continuous and decreasing function taking value in [—1, oo]. 
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p, the resulting SNE is not unique as shown in the following lemma. For notational simplicity we write 

Pt = (eAC) 2 . 

Lemma 4: Assume AC > 1; then there exists a p_i > p t such that for all p satisfying p_i > p > p t 
two SNE are possible. One resulting from the W and the other from the W_i branch of the Lambert 
function. 

Note that for for values of p E [p t , the equilibrium MAP computed on the W_i is larger than 

that computed on the Wo branch. If nodes reach the equilibrium that is resulting from the W_i branch, 
then they will be transmitting aggressively at equilibrium. Later we will see that this leads to inefficiency 
©■ 

Lemma 5: Assume AC < 1, there exists a p > p t such that for all p satisfying p$ > p > pt, the 
SNE is p* — 1. For all p satisfying p > p the SNE is unique and lies on the Wo branch of the Lambert 
function. 

Lemma [3] together with |4] and [5] completely characterize all possible equilibria. We summarize all the 
SNE resutls in Table IB Figures Q] and [2] illustrate the Lambert function and its properties used in the proof 
of Lemma |4] and [5] In Figure \T\ the y axis marked with a double arrow denotes the region in which two 
equilibrium points occur. 

Stability of SNE: In Figure [3] we plot the best response of the tagged node and that of all nodes against 
the tagged node. This example demonstrates the existence of two SNE. The smaller SNE among the 
two results from the principal branch of the Lambert function in (fl9l ), and the larger one from the W_i 
branch. When two equilibrium points exist, we denote the SNE resulting from the W branch of the 
Lambert function as pi and that resulting from the W_i branch as p* x . When we do not need to make 
this distinction or it is unique we write the SNE as p*. From Figure [3] we see that, at the equilibrium 
point Pq, a slight increase in the MAP p results in a decrease in p'. This is a stabilizing behavior and we 
conclude that pi is stable. In contrast, at the equilibrium point p*_ x a slight increase in p is seen to cause 
an increase in the value of p' . Thus the second equilibrium point is unstable. 

Optimal Pricing: Assume that p* is an equilibrium point for a given value of p that satisfies the conditions 
in Lemma |4] or [5] for a given value of A and C. The potential delay experienced by the packets of a tagged 
node at equilibrium can be computed as 
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^ p* exp{— p*XC} 



' (20) 



p* exp{— p*XC} exp{— p*XC} 

Vp_ 

exp{W(-AC/yp)} 

_VpW (-Ag/yp) 
W (-XC/y/p) exp {W (-XC/y/p) } 

y/PW(-XC/y/p) P W(-XC/y/p) 



(21) 

(22) 
(23) 



(-AC/VP) AC 
In above chain of equalities the relation C = 2C is used in (|20l ). (I2TI) follows from (TT71) . Equation (1231 
follows by applying the definition of the Lambert function to the denominator term in (l22l) and rearranging. 
From [23] we see that when the SNE is not unique when the Lambert function takes two values for some 
price factors. Delay experienced by the tagged node is larger when the SNE results from the W_i branch 
of the Lambert function. 

With the expression for the delay of a tagged node at equilibrium, one can look for the value of the price 
factor that minimizes the delay experienced by each node at equilibrium. We assume that the objective of 
the central agent is to minimize the average delay per unit area, i.e., spatial delay density, at equilibrium. 
The optimization problem of the central agent is: 

minimize - pW (z^Z^) subject to yfp > XCe and - W (-XC/y/p) < AC. (24) 

The first constraint in this optimization problem results from Lemma |3] and the second constraint is 
required to ensure that the resulting value of p in (fTTl) lies in the interval [0 1] . 

Let h(p) := — pW i—XC/yfp) denote the objective function in the above optimization problem without 
the multiplicative factor C. h(p) is defined for yfp > eXC. In the following lemma we state some of its 
properties. 

Lemma 6: On the principal branch Wo, h(p) is a quasi convex function in p and the global minimum 
is attained at p* = 4e(AC) 2 . On the W_i branch h(p) is a monotonically increasing function taking value 

(eAC) 2 at p = (eAC) 2 . 

Consider the optimization problem in ([24]) on the Wo branch of the Lambert function. The value of p* at 
which h(p) achieves minima satisfies the first condition in (12"4l) as p* = 4e(AC) 2 > e 2 (AC) 2 . The value of 
p in equation (fTTb at p = p* is 1/(2AC). This implies that whenever 2AC = AC > 1, the resulting value 
of p lies inside (0 1), thus satisfying the second condition of the optimization problem. Hence under the 
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assumption AC > 1, the global minimizer of h(p) lies in the constraint set of f[2~41) and the value of the 
objective function at this point is 

p*W(-AC/v^) (4e(AC) 2 )(-l/2) 



eAC. (25) 
AC AC 

When AC < 1, any value of p such that p < p* violates the second condition, as — W(— XC/^/p) > 1/2. 
From Lemma [6l the minimum value of h(p) is achieved at a p satisfying — W(— XC/yfp) = AC. This 
implies that under the condition AC < 1 the value of the SNE is p* = 1 and the delay experienced by 
each node is given by exp{— AC}. We summarize these observations in the following result. 

Theorem 2: The value of the price factor p* that minimizes the delay at equilibrium is as follows: 

. / (2ACVi) 2 if AO 1 

P = < _ (26) 

[ -2W(-AC/Vp*) = AC, if AC < 1 

and the corresponding delay at equilibrium is 

f A 2 eC if AC > 1 

dt(p\p*)={ (27) 
[ Aexp{AC}, if AC < 1. 

Comparing (1271) and © we conclude the following result: 

Proposition 4: The spatial delay density in the game problem at equilibrium equals the global optimal 
spatial delay density, i.e., d t (p m ) = d t (p*,p*), if the price factor is chosen as in (|26l ). 

Again, by appropriately pricing the nodes, the selfish behavior can be used to attain the global optimal 
performance in the game problem. If the SNE is not unique at the optimal price factor p*, then one 
needs to ensure that the nodes reach the equilibrium that is computed on the principal branch of the 
Lambert function to realize the global optimal performance at equilibrium. Indeed, if AC > 1 and 
— W_i(— l/2y/e) < AC or equivalently p_i > 4e(AC) 2 , at the optimal price p*, an equilibrium point 
exists on the W_i branch of the Lambert function. It would be interesting to learn about the way to reach 
a favorable equilibrium at the optimal price factor. However we do not pursue this question in this paper. 

Remark 1: If AC > 1 and optimization is restricted to the W_i branch in (|24l) . then by Lemma |6l the 
objective function is minimized by choosing a p satisfying Jp = eAC, which results in the equilibrium 
probability 1/(AC). The spatial density of delay at this equilibrium point is given by 

_ P w(-xc /V r) = = 

AC AC 

Comparing this value with ©, we see that spatial density of delay increased by a factor of e/2 by the 
selfish behavior of the nodes. 
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VII. PRICE OF ANARCHY 
In this section we study the degradation in the network performance due to a selfish behavior of the 

nodes. The Price of Anarchy (PoA) compares the social utility at the worst equilibrium with the optimal 

social utility [18]. For our Poisson bipolar MANET with infinitely many players, we define the PoA as the 

ratio of the optimal spatial average performance that can be achieved, to the spatial average performance at 

the worst SNE. Recall that we denoted the system utility by U(p) when we considered the team problem, 

with each node using the same MAP p. For the game problem we denoted the utility of a tagged node by 

U(p',p). In the game problem the spatial average performance at equilibrium is evaluated by multiplying 

the utility of the tagged node and the intensity of the Pp.p . Then the PoA is given by 

max cf =r ii U(p) 

PoA = — pg[ ° 1J r , W - , (28) 

min p . 6 s XU yp*,p ) 

where S C [0 1] denotes the set of symmetric Nash equilibria. 

We study the PoA as a function of p for a given value of A and C. The utilities studied in Section [V] 
and [VI] are considered below. 

A. Goodput 

In this subsection we consider the utility defined in Section [V] Let us begin by considering the team 
utility. When all the nodes use the MAP p, then from Equation (flOl) . the team utility is given by 

U(p) = Xpexp{— pXC} — Xpp. (29) 

If p > 1, the maximum value of the utility is zero and the maximum is attained at p m = 0. Let p m := p m (p) 
denote the MAP that maximizes the team utility in (|29l) . The following lemma gives its value. 
Lemma 7: The MAP value that maximizes the team utility (|29|) is given by 

1 - W(pe) 
Pm ~ AC 

for all p > if AC > 1, and if AC < 1 it is the maximizer for p such that W(pe) > 1 — AC. Further the 
maximum team utility is given by U(p m ) = — " 



By using the definition of the Lambert function, one can verify that U (p rn ) is a decreasing function in 
p. Indeed, differentiating U(p m ) with respect to p we have 



dU( Pm ) (l-W(pe)) 2 (l-W(pe))(l + W(pe))WV) 

pe r — — — r (30) 



dp CW(pe) r CW(pe) CW(pe) 

(1 - W(pe)) 2 _ (1 - W(pe)) 
CW(pe) CW(pe) ' 



(3D 
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In (1301 W denotes the derivative of the Lambert function. Equation (|3T| ) follows by applying formula 
for the derivative of Lambert function. The last equation is negative valued for all p G [0, 1]. Thus the 
optimal utility is a decreasing function in p. 

Let us look at the utility of the tagged node at equilibrium. From Proposition [3] we have 

TT( * J exp{-AC} -p, if p < exp{-AC} 

U{p ,p)=< (32) 
[ 0, if p > exp{-AC} 

The utility at equilibrium is also a decreasing function in p for all p < exp{— AC}. With the expression 

for utility at equilibrium and global optimum we have the following result for PoA: 

Theorem 3: The value of the PoA is 



p(l-W(pe)) 2 



if p < exp{-AC} 



PoA(p)= { ACW(p e){e xp { -AC}-p}> - - - (33) 

oo, if p > exp{— AC} 

The PoA is shown as a function of p in Figure [5J From this figure we see that as p increases, the PoA 
grows unboundedly. Thus the PoA is optimal when the pricing factor is set to zero. If AC > 1 then the 
PoA is infinite by definition at p = 1/e. However we noted in Section [V] that the optimal performance of 
the spatial density of success is achieved at equilibrium with the same price factor. If AC < 1 then the 
PoA is infinite at p = exp{— AC}. But again we noted in Section |V] that at this price factor the optimal 
performance of the spatial density of success is achieved at equilibrium. 

In Figure[6]the equilibrium MAP and global optimal MAP are shown. For all values of p, the equilibrium 
MAP is larger than the global optimal MAP. Hence the nodes transmit more aggressively at equilibrium. 
But we note from Figure |6] that the gap between the global optimal MAP and the equilibrium MAP 
reduces with pricing. 

B. Delay 

Consider the utility function in Equation ([Tot . The team utility for this game, when each node transmits 
with MAP -p is 

U( P ) = ^V7n ~ VP- (34) 

It is easy to verify that the above utility function is concave in p. Assume that AC > 1. Then the unique 
MAP, denoted as p m := p m (p) that maximizes the social utility satisfies 

exp{p m AC}(l - p m AC) = pp 2 . (35) 

We obtain this by differentiating Equation (|34l ) and setting to zero. Note that any p m that satisfies Equation 
(1351) also satisfies p m \C < 1, hence p m € [0 1]. Also, it can be easily verified that p m is decreasing in p. 
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The utility at equilibrium can be obtained by using the equilibrium MAP in ([191 ) and (1231) , as 

U(p*,p*) = (-2\p/\C)W (^] ■ (36) 



V VP 

When AC > 1, from Lemma @] two symmetric Nash equilibria are possible for the price factor p < p_\. 
Hence the above utility function can take two values, one corresponding to each equilibrium. Recall that 
we denoted by p^ the SNE computed on the principal branch, and by p*_ 1 that computed on the other 
branch of the Lambert function. Recall that pi < p*_ v The following proposition gives a bound for the 
PoA 

Theorem 4: For Poisson bipolar MANETs with utility as in(fT6l). 



Pm(p-i)A(2 -pm(p-i)XC) n p m {p t )\C{2 - p m {p t )\C) 

2(l-p w (p_ 1 )AC) " ° {P) - 2(1 -p m ( Pt )XC) K ) 

for p E [p t p_i], where y/p^ = eXC. In addition 

^ < PoA(p) < 1 for p > p_ x . (38) 
Po 

The PoA as a function of p and the bounds obtained in Proposition |4] are shown in Figure [7J The jump 
in the figure at p = p_i is due to two possible SNE for p < p_ x and a unique SNE for p > p_ x . In the 
interval [p t , p_i] the PoA is decreasing in p. This results from the bad Nash equilibrium that occurs on 
the W_i branch of the Lambert function which increases in p. If the central agent can't set a price factor 
higher than p_i then, from the PoA point of view, it is desirable to set the lowest possible price factor, 
i.e., p = Pt- For a price factor larger than p_ 1? there is unique SNE, which is smaller than the equilibrium 
that occurs on the W_i branch and decreases^ with p. Thus setting high a price factor leads to improved 
PoA. 

If p_i < 4e(AC) 2 , from Proposition 0] by setting p = 4e(AC) 2 , we can obtain better performance 
at equilibrium and also good PoA. If p_i > 4e(AC) 2 then by setting p = 4e(AC) 2 one obtains global 
optimal performance at equilibrium provided the nodes settle at an equilibrium that lies on the principal 
branch of the Lambert function, otherwise this price factor leads to a poor PoA. 

VIII. Conclusions 

Geometric considerations play a very central role in wireless communications, since the attenuation of 
wireless channels strongly depend on the distance between transmitter and receiver. Models that take into 
account the exact location of mobiles are often too complex to analyze or to optimize. Our objective in 

'Principle branch of Lambert function is decreasing function of p 
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this paper is to model competition between mobiles as a game in which the locations of players is given 
by a Poisson point process. 

More structured point processes can also be contemplated, for instance exhibiting attraction (hot spots) 
or repulsion (more elaborate medium access control than Aloha like e.g. CSMA). We leave the analysis 
of medium access games under such point processes for future research. 

The competition we considered in the paper was between individual mobiles each taking its own selfish 
decisions. We saw that the equilibrium of the game results in a more aggressive access (larger access 
probabilistically). We studied further pricing, and identified pricing parameters that induce an equilibrium 
achieving the social optimal performance. On the other hand we showed that the utility at equilibrium 
can be zero). 

We plan in the future to study other games within this framework: for instance games with finitely 
many operators, each taking decisions for all its subscribers. In addition we shall study jamming games. 
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APPENDICES 

Proof of Lemma |4] 

First consider the Wo branch of the Lambert function. As p takes value in the interval [(ACe) 2 oo], 
W increases continuously from —1 to 0. Thus (TT71) has a solution in the interval [0 1] if AC > 1. 

This implies that equilibrium point exists on the W branch for all p, satisfying yfp > ACe. 
The W_i C~jjf) branch decreases continuously from — 1 to — oo as p takes value in the interval 
[(ACe) 2 oo]. This implies that there exists a sfp_ x > ACe such that W_i C==^-j) = -AC, and for 
all p such that y / p_ 1 > y/p > Ce satisfies — W_i C~jjf) — resu l tm § in a p G [0 1] that is a solution 
of (fTTT) . Hence there exists an equilibrium point on the W_i branch for all p satisfying • v /p_ 1 > y/p- This 
concludes the proof. 

Proof of Lemma \5\ 

As in Lemma @] we can argue that on the Wo branch, there exists p > ACe such that Wo ( — 

V V Pa 

— AC and for all p such that p > p_ x > ACe satisfies — W \ — XC as — W is decreasing in p. 

Proof of Lemma |6| 

Differentiating h(p) with respect to p 

±h{p) = -W {-XC/y/pj - W (-AC/Vp) (AC/2VP) (39) 

= -W (-XC/y/p) | 1 r- I . (40) 

1 /V; \ 2(1 + W(-AC/Vp)); 

In Equation (l39l) W'(-) denotes the derivative of the Lambert function which is given as ifTTl reqn. 3.2] 

W ^ = nYwr for x * °' x * ~ 1/e - (41) 
x{l + W(x)) 
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Equation (|40l ) is obtained by substituting the derivative in (14T1) . evaluated at x = —XC/y/p, in Equation 
(|39| ). Recall that on the principal branch of the Lambert function W(— XC / yfp) is a negative valued 
increasing function in p. Then the term within parenthesis in (|4Q|) is a increasing function of p passing 
through the origin at p* that satisfies W(— XC / y/p*) = —1/2. Thus h(p) is decreasing for p < p* and 
increasing for p > p*. From [19J[sec. 3.4.2] we conclude that h(p) is a quasi convex function in p. 
Further by the definition of the Lambert function 

-XCjyf^ = W(-AC/v^)exp{W(-AC/v^)} = --exp{-l/2}. (42) 

2 

Rearranging Equation (|42l) . we get p* = 4e(AC) 2 . The other part of the Lemma follows by noting that 
— W(— XC/y/p) is an increasing function in p on the W_i branch. 

Proof of Proposition M 

From Equation (13~6l) and (|34|) we have 



(44) 
(45) 



(-2p/AC)W (-^ 
p m /(l -p m AC) +p r , 



;-2/AC)W(-^j 

p m (2 -p m AC) 



(-2/AC)w(-^) (l-p m AC) 
(-l/AC) W (-f) 



> Pm (46) 



We arrive at equality (1441) by dividing both numerator and denominator in (1431 by p, and applying the 
relation in (T351 ). Equality (1431 ) is obtained by simple rearrangement of terms in the previous step. 
To derive the bounds in (1371 ), we consider the equilibrium computed on the W_i branch of the Lambert 
function as it leads to the worst case equilibrium. This equilibrium is an increasing function in p in the 
interval po < p < P-i as discussed in the proof of Lemma @] Also, recall that the value of p m is decreasing 
in p. Thus the numerator in (1441) is decreasing in p. Which implies that the ratio in (1431 ) is also decreasing 
in p. The upper bound in (I3~7l) now follows by noting that — W(po) = 1- To obtain the lower bound we 
use the relation -W(-AC ? / v /pTT) = XC in (I4~6l 

For values of p larger than p_\ the SNE is unique, resulting from the principal branch of the Lambert 
function. The upper bound in (|3~8l) follows directly by the definition of PoA, and the lower bound follows 
from the Inequality (|46|) . Note that the lower bound is a function of p. 
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Proof of Lemma [2] 
Consider a point to set map 5 : [0, 1] — > [0, 1] defined by 

<J(p) = ^p'|C/(p',p) = max U(q,p) > . 

This defines the set of best responses of the tagged node when all the other nodes use the MAP p. It 
follows from the continuity of U(p',p) in p and concavity in p' for a fixed p, that T is an upper continuous 
mapping that maps each point of the set [0, 1] into a subset of [0, 1]. By the Kakutani fixed point theorem, 
there exists a point p* e [0, 1] such that 

U(p*,p*) = max U(p',p*). 
p'e[o, l] 

Then p* is the symmetric Nash equilibria by Definition [TJ 
Proof of Lemma [7| 

Assume p < 1, then p m > and satisfies the relation exp{— p m AC}(l — p m \C) = p. By rearranging, 
p m can be expressed p m = (1 — W(pe))(AC). Recall that W(-) is a monotonically increasing function 
taking values W(0) = and W(e) = 1. If AC > 1 the p m lies in the interval [0, 1] for all p e [0, 1]. If 
AC < 1, then p m lies in the interval [0, 1] for all p such that W(pe) > 1 — CA. Thus whenever AC < 1 
we assume that p satisfies W(pe) > 1 — CA. The optimal value of the utility function is 



p m exp{p m AC} = 1 ^ (Pe) exp{-l + W(pe)} (47) 

' {exp{W(pe)} -pe} (48) 



eAC 

1 I " ,„1 "I 1 , 



eAC\w(pe) p6 j AC \W(pe) 1 j (49) 
where (1471) is obtained by substituting the value of MAP maximizes the team utility. Equation (1481) and 
(l49l) follows by application of the definition of the Lambert function. The maximum utility for the team 
case can be now computed as a function of p 

p(l-W(pe)) 2 

U{Pm) = CW(pe) • (50) 
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Table I 

Characterization of symmetric Nash equilibria 



AC 


Price factor 


SNE (p*) 






P < Po 


1 




AC < 1 


P > Po 


-W (-AC/^p)/AC 






p < (eAC) 2 


1 




AC > 1 


P-i > P > (eAC) 2 
P > P-i 


-W (-AC/^/p)/AC or -W_i(- 

-w (-ac/vp)/ac 


\C/^p)/\C 
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Figure 1. SNE on both branches 
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Figure 2. Unique SNE 




Figure 5. PoA for Goodput 



Figure 6. Equilibrium and Optimal MAP 




Figure 7. PoA for delay, C=3, A=l 



